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We study the dynamical effect of strange quark mass as well as kincmatical one on 
the color- flavor unlocking transition using a NJL model. Paying a special attention 
to the multiplicity of gap parameters, we derive an exact formula of the effective 
potential for 5-gap parameters. Based on this, we discuss that the unlocking 
transition might be of second order rather than of first order as is predicted by a 
simple kinematical criterion for the unlocking. 



1. Introduction 

The dynamics which quark-gluon matter exhibits under high baryon den- 
sity is one of the most challenging and exciting problems in QCD 2 ' 3 - 4 . A 
number of literatures have revealed the realization of the color-flavor locked 
(CFL) type of the pairing order 5 for sufficiently large quark number chem- 
ical potential 4 > 6 < 7 . 8 < 9 > 10 . I n contrast to this solid fact, however, the phases 
for large, and realistic value of M^/fi are still veiled in mystery, and so 
many phases and the associated phase transitions have been studied so far 
Ii,i2,i3,i4,i5,i6,u Thc unlocking transition from the CFL phase to the 2SC 
state is simplest one of the examples. This transition with increasing the 
value of M s 2 /fj, is firstly investigated using Nambu-Jona-Lasinio (NJL) type 
effective models 11:12 . Their results show that the 1st order unlocking tran- 
sition takes place at some critical , and a simple kinematical picture for 
the critical mass works in quite satisfactory way. This criterion is based on 
the conjecture that the transition occurs when the mismatch of the Fermi 
momenta of light and strange quark becomes as large as the magnitude of 
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the gap. In Ref. 11, the expansion in flavor breaking parameter M^/jj, is 
used to the construction of the fermionic dispersions and the effective po- 
tential. On the other hand, in Ref. 12, the coupled gap equation is exactly 
solved with non-perturbative treatment of M s . However, their construction 
of the effective potential is quite ambiguous and, as we shall show, possibly 
fails because of the multiplicity of gap parameters. Even though their re- 
sults are in good accordance with each other, there still exist a possibility 
that the true solution is missed. 

In this talk, we revisit the unlocking transition in more complete way 
than others 11 ' 12 1 by making a proper use of the Pauli construction of the 
effective potential. In particular, we study how the CFL state and other 
states are realized in the multi-gap parameter space and how the potential 
gets distorted with increasing the value of M| / fi. More complete analyses 
have been made in Ref. 1. 

2. Coupled Gap Equation and Effective Potential 

In this section, we present an outline for obtaining the gap equations and 
the effective potential with which we can determine the ground state of the 
system characterized by (/j,,T,M s ) parameter. 

Self-energies. We first introduce the color-flavor mixed quark base: q a = 
J2 a i (^"/ V^Ya^i' wnere A Q (a — 1 <~ 8) are the Gell-Mann matrices and 
we defined (A 9 )" = yJ^Sf. Then the Nambu-Gor'kov propagator in the 
2-compnent quark field Q = (q, (f) is written as 



S(q ,q) = iF.T.(T{Q(x)Q(0)}) = 




T.(qo,q) is off-diagonal self-energy which gives rise to gaps in the quasi- 
quark dispersion. M is the quark mass matrix in the color-flavor mixed base 
M a p — tr [\ a rh\p] /2, where m = diag.(0, 0, M s ) is mass matrix in flavor 
space. The pure CFL ansatz for the off-diagonal self-energy is expressed in 
the color-flavor mixed base as following, 

^o, 9 ) = C 75 0( A8l8 Ai ). (2) 

Here, C is the charge conjugation matrix C — 17270- C75 guarantees 
that the pairing takes place in the J p = + channel. On the other hand, 
we would have the 2SC state with u-d pairing £ 2SC " d = -e i:) e afc3 A2sc = 
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AA^ije d for sufficiently large value of M s , which is written as 



/A 



2SCJ-3 



\ 



o 4 



(3) 



\ -— A 2 sc -3A2SC j 

in the color-flavor mixed base. From the two expressions for the 2SC and 
CFL phases, it is quite natural to expect the distorted CFL state (dCFL), 
the minimal interpolating pairing ansatz between those phases, for the small 
but finite strange quark mass 11 ' 12 , 



v dCFL/„ \ 



M s <M r : 



C75 



/A 83 ] 



V 



A82I4 



\ 

A 8 i A 
A A n / 



(4) 



In this parameterization, the SU(3)c+v symmetric CFL state and the 2SC 
phase with SU(2)c x SU(2)l x SU(2)r symmetry are expressed only as 
different limits of the g?CFL phase which is invariant under color flavor 
simultaneous rotation SU(2)c+v- By introducing the 5-dimensional vector 
A = (A 83 , A 8 2, A 8 i, A, An)*, we can express the model space for the pure 
CFL state as Acfl = (A 8 , A 8 , A 8 , 0, Ai)*, which spans a 2-dimensional 
planer section, while the 2SC phase as the I-parameter vector A 2 sc = 
(A 2SC , 0, -A 2SC /3, -y2A 2SC /3, "2A2SC/3) 4 . 

Gap equations. The anomalous propagator is defined by the off-diagonal 
element of the Nambu-Gor'kov propagator, which takes the following form 
for our ansatz Eq. (4), 



S^(q)=iF.T.(q a (x)qe(0))= 75 C< 



R82I4 



V 



\ 

Rsi R 
R R\\l 



(5) 



{R&3, R&2, R&i, R, Ru} are propagatores, which are complicated functions 
of (q, A, /i, M s ) 1 . The gap equation is obtained by using Feynman rule 
for the self consistency between proposed self-energy and the one-loop self- 
energy. In the case of NJL model with interaction vertex for one-gluon 
exchange, £; nt = (l/2)g 2 q'f fl t a qq-f^t a q 7 we obtain 



Sa/3 = iVy ^(T% a S$(q)(T a ) 50 . 



(6) 



February 1, 2008 23:8 WSPC/Trim Size: 9in x 6in for Proceedings 



H.Abuki 



4 



Here, the bare vertex in the color-flavor mixed base is defined by (T a ) Q( g = 
tr [AaAaA/j] /4. This integral matrix equation contains the following set of 
five equations. 

A 83 = (57e 83 - K 81 - 2^2K - 272-n) = g 2 K &3 [A], (7) 
A 82 - g 2 -^ (272. 8 2 + 272: 8 i + V272- - 272-n) ee g 2 K S2 [A], (8) 
A 81 = g 2 -^ (-37e 83 + 872. 82 - 72: 81 + 2V272. - 272-n) = g 2 K sl [A], (9) 

A = g 2 ^- (-372. 8 3 + 272. 82 + n S i - v/272.) ee g 2 K[A], (10) 

An =g 2 \ (-37l 8 3- 4-R 8 2 -7e 81 ) = g 2 K n [A}. (11) 
Here, we have defined 72- by 

72-i(A; M s ) = i J Jtr [c 75 R t ( q] A, M s )] . (12) 

Effective potential for multi-gap parameters. Here, we do not attempt to 
formulate the Pauli-construction of the effective potential, which is exactly 
done in Ref. 1. But instead, we illustrate how we would miss the true 
effective potential due to the multiplicity of gap parameters when it is 
constructed by the integration of the gap equation, and only show the 
correct procedure to obtain true one. 

We might think that the derivative of the effective potential is 

^° ff oc (gap equation for Ai). (13) 

If so, the effective potential is constructed by 



(14) 



But this gives a false formula in the case that the many gap parameters 
exist and couple each other by different couplings. 

Actually, the derivative of the effective potential coincides with some 
linear combination of the gap equation like 

9 2 DT j ^=A ij (A j -g 2 K j [A}), (15) 

where the diagonal matrix D = diag.(3, 4, 1, 2, 1) represents the degen- 
eracies of the gap parameters, and the matrix A is relating the gap pa- 
rameters A = (A 8 3, A 8 2, Agi, A, An) and the corresponding condensates 
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= (083: 082, 081, 07 0n) as 4g 2 = AA. For the detail of the definition <j> 
and the form of the matrix A, see Ref. 1. We can obtain the true formula 
for the effective potential as 

n cS [A] = ^^AA- -£ / dt{DA) t3 A t K 3 [tA\. (16) 
DA is real-valued symmetric matrix, whose eigenvalues are all real. 
3. Numerical Results 

Here, we present our numerical results. The cutoff parameter is set to 
A = 800MeV, and the coupling parameter g 2 is tuned to reproduce 400MeV 
for the constituent mass of quark at zero chemical potential 12 . 

3.1. Solution of gap equation 

We display the solution of the coupled gap equation Eqs. (7)~(11). In 
FIG. 1(a), we show the eigenvalues of the original gap matrix Eq. (4) as 
functions of M s . (A 83 , A 82 ) are the gaps for iso-triplet and doublet modes. 
(XiiX2) are the eigenvalues of the iso-singlet mixing sector. According to 
the values of these parameters, the states are distinguished by 



Gap parameters 


Global Symmetry 


States 


A 83 = A 82 = X2, Xi 7^ 
A 83 = Xi, A 82 = X2 = 
otherwise 


SU(3) V 
SU(2) L x SU(2) R 
SU(2) V 


CFL for M s = 
2SC for M s > M c s 
distorted CFL (dCFL) 



Continuous Color-Flavor Unlocking. It is surprising that the system seems 
to undergo continuous transition from the CFL to the 2SC, and actually 
stays in the distorted CFL (dCFL) phase even for M s > a = 400 MeV. 
We will see later that even in our case, the 2SC is always a solution of the 
gap equation, but with higher energy than the CFL state. In FIG. 1(b), 
we display in FIG. 1, the gaps (A 83 , A 82 , Xi> X2) for M s — 250 MeV as 
functions of temperature. We see that as the temperature is raised, the 
dCFL phase encounters the continuous transition to the 2SC state at T = 
35MeV, and eventually the system undergoes 2nd order phase transition to 
the normal fermi-liquid at T — 60MeV. 

3.2. Effective potential for multi-gap parameter space 

We now study how these states are realized in the multi-gap parameter 
space by computing the effective potential with help of the formula Eq. (16). 
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(a) (b) 




100 200 300 400 500 600 10 20 30 40 50 60 70 



MJMeVJ r[MeV] 

Figure 1. (a) Gaps are plotted as functions of M s . (b) The temperature dependence of 
the gaps for M a = 250 MeV. 



Let us first introduce the vector 

_ (,-10)^4-0-10)^ 
50 

with the CFL gap vector A CFL = (80,80,80,80,-175) and the 2SC gap 
vector A 2S c = (106,0,-35,-50,-70) for Ms = (in the chiral limit). 
Eq. (17) define 2-dimensional planer section in the 5-dimensional gap pa- 
rameter space, which includes the simple Fermi gas (10, 10) with A = 0, 
the CFL (60, 10) and the 2SC (10, 60). In FIG. 2, we show the contour plot 
of the effective potential Q qS [A(i, j); M s ], for M s = MeV (a), 200 MeV 
(b), 300 MeV (c) and 400 MeV (d). We can see that the CFL state (60, 10) 
looks the global minimum in this 2-parameter space (i, j) for M s = MeV. 
On the other hand, the 2SC state (10, 60) is realized as a saddle point which 
is stable in the direction of the simple Fermi gas (10, 10), but is unstable 
in the direction towards the CFL state (60, 10). As the strange quark mass 
is increased to 200 MeV, the CFL minimum moves towards the 2SC point, 
and its condensation energy gets reduced, while the position and the en- 
ergy of the 2SC state is unaffected by M s . This minimum point (55, 20) is 
expected to be located close to the distorted CFL dCFL state in the full 
5-parameter space. The dCFL state gets distorted significantly towards the 
2SC state at M s = 300 MeV (FIG. 2(c)), and seems absorbed into the 2SC 
as the strange quark mass approaches the order of the chemical potential 
~ [L of the system. We can conclude that the unlocking transition is not of 
1st order. More detailed analyses made in Ref. 1 reveals that this transition 
is of 2nd order and indicates that the g?CFL satate at low strange quark 
density is the Bose-Einstein condensation of tightly bound pairs rather than 
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(a)Ms=0[MeV] (b) Ms=200 [MeV] — 1 




[MeV/fm 3 ] 



Figure 2. Contour plot of the effective potential Q c ff [A(i, j)] where A(i, j) is defined 
by Eq. (17), for M s = MeV (a), M a = 200 MeV (b), M s = 300 MeV (c), and for 
M s = 400 MeV (d). 

the BCS state 18 . 

2SC as a solution of gap equation. What should be stressed here is that 
the 2SC stays at least a saddle point in the full model space irrespective of 
the value of M s . This means that the 2SC is always a solution of the gap 
equation, and if we had missed in obtaining the correct effective potential, 
then we misunderstood it as the true ground state instead of the dCFL 
state for intermediate strange quark mass M s . 

Unstable CFL state? In FIG. 3(a), we plot the section of the effective 
potential surface FIG. 2(a) on the line linking the 2SC (10, 60) A t= o = ^2SC 
and the CFL (60, 10) A t= i = Acfl- Both states are determined by solving 
the gap equation in 5-gap parameter space for M s = MeV. Now we address 
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Figure 3. (a) The values of the effective potential on the 1-dimensional line con- 
necting from the 2SC state t = to the CFL state t = 1. (b) The bold line is 
just an enlargement of the potential curve for M a = in (a), and five set of the 
points (triangles, bold dots, squares, reflected triangles, lozenges) represent the po- 
tential curves in the directions of the eigenvectors belonging to the five eigenvalues 
A ^ (0.19, 0.14, -0.99, -1.70, -2.20)/* 2 of the Hessian curvature matrix at the CFL. 



the question whether the CFL state is true global minimum even in the full 
model space or not. FIG. 3(b) shows us the answer for this question. Five 
potential curves are shown near the pure CFL for M s = 0, each of which 
represents the potential curve in the direction of an eigenvector belonging 
to the corresponding eigenvalue A of the Hessian curvature matrix at the 
CFL point d £.Q A _ _ ■ What is surprising is that the CFL state is not 

1 j A=Acfl 

the global minimum in the full model space. It is unstable in three out of 
five directions. These three directions correspond to the color-flavor sextet 
channels in which the interaction acts repulsive 1 . Also it should be noted 
that the color-flavor sextet gap parameters do not break any symmetry, 
thus are not the order parameter for the unlocking as is anticipated in 
their critical behaviour ~ (1 — T/T c ) 3 / 2 near critical temperature T c in the 
chiral limit 19 . Anyway, because the effective potential becomes unbound 
once the sextet condensations are included, it makes no physical sense to 
include symmetric mean fields contribution 5 > 20 > 21 in the ansatz for the gap 
matrix. It is said that the sextet condensations are not induced within the 
mean field approximation, whereas they might be triggered by higher order 
fluctuation effects beyond the mean fields. 



4. Conclusion 



We have adopted the NJL model to study the phases of quark matter under 
high chemical potential. By making proper use of the Pauli-construction 
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method, we have derived the exact formula for the effective potential for 
multi-gap parameters. In particular, we have studied the unlocking phase 
transition from the CFL to the 2SC. We list main results below. 

1. Second order phase transition. The unlocking transition is of con- 
tinuous weak 2nd order. The 1st order unlocking never appear even at 
T = 0. This contradicts the simple kinematical picture for the color-flavor 
unlocking transition. 

2. Toughness of the CFL state. The CFL state at T — is much more 
robust against the increase of M s than is predicted from the kinematical 
criterion, and the 2SC state is continuously connected from the dCFL state 
at the strange quark mass > /i 2y/ /iA§(/i, M s = 0)). 

3. The 2SC as a saddle point. 2SC state is always a saddle point, a 
solution of the gap equation, of the effective potential for any value of 
M a , which is unstable in the direction towards the dCFL state, as well as 
three color-flavor sextet directions. The (iCFL state, a solution of the gap 
equation with a larger condensation energy than the 2SC state approaches 
the 2SC saddle point as M s approaches M£ . 

4. Role of the sextet gap parameters. The dCFL state is also unstable in 
the directions for three color-flavor symmetric channels. This is attributed 
to the fact that the interaction is repulsive in these channels. We should 
not include the symmetric components of the gap parameter into our ansatz 
from the beginning, because no Cooper instability is present in these chan- 
nels due to the absence of the attractive force. However, the effects of 
the sextet gaps on the anti-triplet sector are very small even if they are 
included. 

In this talk, we completely neglected the electric charge neutrality and also 
the color neutrality 22 . It would be interesting to include these effect into 
the gap equation, and to study how our dCFL phase is robust or fragile 
to be withdrawn by the neutrality condition. Also, we have ignored the 
usual chiral condensates in the vacuum. However, we have to include this 
effect to discuss the phase transitions in the lower chemical potential region. 
Studying phase transitions in this regime by improving our model along this 
line is also an interesting subject to be done in the future. 
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